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As a general rule the following conventions will be observed in 
this report. A superscript (T) will denote a total quantity, i. e. the 
result of superposing initial pressurization effects and other loading 
effects. A superscript (n) will denote the r» ia verm of a Fourier ex- 
pansion of the quantity. (See Art 3. U). A "primed" quantity will de- 
note a nondimensional quantity. There are two exceptions to this last 
rule, E' and X . For explanation of these exceptions and the method 
of nondimensionalization, see Art. 4. d Certain terms will carry a 
double sign, e. g., +■ , the "upper sign" should be used if the "upper 

series" expansion of Art. 3. 0 is usee, 

The symbols have the following me am a js: 

<t> , 0, z - - -coordinates of a point of tne shell; z is measured 
positive inwar d from the middle surface, q and 0 
are the usual angles of a spherical coordinate system, 
u , v, w components of the displacement 

^4><t>* ^00' ^ <!>(/ %z‘ N oz - - -internal stress resultants 

M gg, - — internal stress couples 

X , Y , Z - - -components of surface traction 

th 

n - - -denotes the n term of tne expansion 

E Young's modulus 

v - - -Poisson's ratio 

a - - -undeformed radius of the middle surface 

h thickness of the shell 

- - -Fourier coefficient of the temperature increase 
• IS - - - Coefficient of thermal expansion. 
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In addition, certain symbols are used which do not appear in the 
hnal equations. These symbols are defined at the time of their intro- 
duction. 
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Chapter 1 . 

1.0 Introduction 

In this report a system of ordinary differential equations govern- 
ing the behavior of spherical shells under arbitrary loads is derived. 

The resulting system is presented in the chapter as equations 1. 1. 1 
through 1. 1. 16. The form of these equations is such that a numerical 
integration technique, e. g. Runge-Kutta. can be easily used. The 
actual method of solution will be discussed in Chapter 5. 

The accuracy of any engineering analyses ts limited by the start- 
ing equations. The basic equations used m this report are felt to be 
the most accurate available at present, l.inear equations have been 
used, i. e. products of displacements, etc. are presumed negligible. 
However, this analysis is designed to include the effects of initial pres- 
surization, and in this sense provides on important extension. In addi- 
tion, the only further assumptions are those of the validity of Hooke's 
law and the common Kirchhoff hypothesis of thin shell theory. 

Several advantages result from this t . pe of analysis. First, the 
basic equations are relatively free of simplifying assumptions. Second, 
very complicated loading systems can be handled. Third, the "output" 
is in the form of the quantities of most use to the designer, i. e. , stress 
and displacement. 

In addition to the assumptions listed above, one additional restric- 
tion has been added. The thickness and material properties have been 
assumed constant. However, it is possible to use the equations 1. 1. 1- 
1. 1. 12 in the special case where the loading is axi symmetric and the 
thickness and material properties vary only in the direction of the gener- 
ator. A slight modification of the equations will allow variation of thick- 
ness and material properties in the direction of the generator for the 
case of a nonsymmetric loading. 
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1.1 First Order F or m_of_ Equation s 

The following eight equations are the i\ nai form of the equations 
which are, for a particular problem*, to he integrated numerically. It 
should be noted that these equations are nondirnensional. (See Chapter 
4 t Section 2, ) The effects of temperature gradients are included in 
these equations. 
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In addition to the eight dependent variables computed above, 
the following quantities are of interest as "output data". 
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Chapter 2. 

2 . Derivat ion of Equations 

The basic equations used in this report are based on the work 
f Dove (A Treatise on the Mathematical Theory of Elasticity, 4th. Ed. , 

: ,*.vr Publications). The rationale of l.ovo’s work will not be discussed; 
.vcV’. r, it may be pointed out that it is felt that his equations are the 
: available. One deviation will be made from Love's work. In the 
• aressions for the stress- displacement relations, we use the work of 
,, (Applied Elasticity, 1st Ed., McGraw-Hill). It was found that 
•,y..r.g retains one more term than Love. The retention of this extra 
term seems justified in that other terms of the same order are present. 
Love's analysis. 

.Since this analysis includes the effect of initial pressurization of 

(T) 

tr.e sphere, the notation N etc. will b usees. The superscript (T) 

Y V 

.:uheates that this is the total stress resultant and can be written as 

N = N -1- N 2. 0. 1 

bb bb 

A'.i'ig N is the stress resultant due to the initial pressurization and 

is the stress resultant caused by loads other than the initial pres- 
s'. . Ali other quantities use a similar notation. 
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2. 1 Equilibrium Equa tions 

On page 535, Love gives the equilibrium equations as: 
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Equations 2.1.1 - 2. 1. 6 are extremely general. We introduce the 
following specialization to a spherical shell: 
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For an explanation of the meanings of the various terms of Equa- 
tions 2. 1. 7, see Figures 2. 1. 1, 2. 1. 2 and 2. 1. 3. 

On using Equations 2. 1. 7 in the equations of Love, p. 523, we find: 
, cos* aw< T > , i a 2 w< T > 

Pl 3 aTsin* ^ aSin * 


2 . 1.8 


p 2 = Sin * + a Sin* 


1 a 2 W (T) , 1 a V 


(T) 


*1 


V = 


, Cos * 

a (T) 
a w 

Cos* (T) 

+ a 

a * 

+ a - u 

1 . i_ . 


1 »u< T > 

a 

a ** 

a a * 

i a 2 w 

(T) x 

a cot * a w 

a a 0 a * a 

TU + a a "6 






2. 1. 9 


2 . 1 . 10 


Cos * (T) _ Sin* 8 v 


(T> 


2 . 1 . 11 




FIG. 2.1.1 


u,v, and w are the components 
at displacements in the direc- 
tions of T,T and"k respectively 
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For a uniformly pressurized sphere with additional loads, we note: 
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On substituting Equations 2. 1. 7 - 2. 1. 15 in 2. 1, 1 - 2. 1. 6, we find 
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2, 2 S t ress- Disp lac ement R el ations 

The stress resultants and stress couples (see, for example, 
Wang, p. 340) are 
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Introducing 
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and using Equations 2. 1. 14 and 2. 2. 2 in 2. 2. 1 we obtain: 
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Chapter 3. 


3. 0 Reduction to Ordinary Differential Equations 

Since the solution, numerical or analytic, of the systems of 
partial differential equations derived in Chapter 2 is very difficult and 
time consuming, it is expedient to expand the dependent variables in a 
Fourier series in the independent variable 9 . This expansion reduces 
the system of partial differential equations to a Bystem of linear ordin- 
ary differential equations which are much more tractable. 

We now postulate expansions of the dependent variables in 
trigonometric series of the form 
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where the symbol { g?^ indicates either a Cosine series or a Sine 

series used. In the following equations certain terms will have a double 
sign, e. g. + ; this notation indicates the "upper sign" should be used 
when the "upper series" is usecl etc. 

Substituting equations 3. 0. 1 in equations 2. 1. 16 - 2. 1. 20 and 
2.2.4* 2. 2. 9 and dividing out common terms, we obtain the following 
system of ordinary differential equations. 
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Chapter 4. 


4. 0 Reduction to First Order Form 


In this chapter, the ordinary differential equations, 3.0.2 through 
!,.(). 'i&, mV iv duced iO ft 't&vth 'pfti’Uctflttt’fy suited to mmitmiud *itw:p’ir- 
tion. The reader primarily interested in the final form of the equations 
is directed to Chapter 1, equations 1.1.1 through 1.1. 16, where the 


1 uc key to successful tmuiipuluUon of these equations lies in the 

introduction of three new dependent variables, X^ n \ <{/ ^ and ^ 

(n ) o 

The quantity \ represents the slope of the median surface in the direc- 

i \ 

tion of the generator. ^ ^ ' is a quantity reminiscent of and related to 

the classical Kirchhoff shear of plate theory. The form of the quantity 
is suggested by the equations and is related to the effective mem- 
brane edge shear for shells of revolution. Throughout this chapter, the 
main purpose will be to manipulate the basic equations into a form such 
that on the left side of the equation we have a single first derivative 
while on the right side we have only the quantities themselves. This 
particular form is conveniently integrated by the Runge-Kutta method 
or some similar integration scheme. 

We introduce the new dependent variable 


x <-» . I 

a d <J> 

Eliminating between 3.2. 3 and 3.2.5, we get 


4.0. 1 



4. 


4.0. 2 


dN ei n) : dM e ^ OCoti {N *“) ij M < n > 

“3* + I d<T"“ + 2Cot *^ + t M e* 


— n *, (n) L n m/r ^ r d? — n 

+ -Emf N 00 - T5TH* M 00 + 1 ’ d ^2 + "Sli 


+ v fc) ] ii + a Y (n) - 

J a 


T “ ...1 

+ ^IHT w 


Eliminating Ng z ' n * between 3. 0. 4 and 3. 0. 5, we obtain 


d ^z (n> - „ dM 6 A (n> 

3$ a5in^ d<p 


+ Co1 * N *z <nl + “l* 1 " 1 + N »+ <n> 


dX to> , N du* n> _ n 2 N ,..(n) . nN „(n) 

■ar- + T Sf- - ^ 


aSin 0 


+ N Cot 0 

x (n) + N_ 
a 

' co. 4 » 

9 ; 

— 

— 2nCoto 
aSin0 

■ M (n) 
0 0 

+ aZ W 


We now introduce a second new dependent variable, 
+ . (n> * asm* K M + i>M, t W 


which implies 

dij> ^ N, *“) < 

~4~ “ aSin * d -Sf— + aCos * N 0z + n 


Substituting 4. 0. 5 in 4. 0. 3 we obtain 




4 . 0.3 


/ 




(ri> 


9_ 


a Sin * d* 


+ N M (n>+N M ln) + N 


(n) 


„<"> + »*L + . W + NCot* A. 


(n) 


C TT, 

a Sin * 


+ 2- cot * u W 


+ a - 0 


n 


a Sin 1 





4.0.6 


Introducing the third new dependent variable 


+ 9 <n) ■ a Sin 4 + Sin* + 


(n) 


4.0.7 


we also get 
(n) 


d.|i 


0 


dN 


d* 


a Sin* 


+ Sin * 


(n) 

— + a Cos0 N 0( ^+Cos* 


~\ 

(n) 


dM 


(n) 


0 * 


4.0.8 


Substituting 4. 0. 7 and 4. 0. 8 in 4. 0. 2 we obtain: 


d4* 


(n) 


0 


T*" 


- {NSin*} v (n) ± {nN} + {na} 


(n) 


2„(n) 


+ {n} M - {Cot*} »J' fl (n) "• { NSin * } — — 


4 00 

- { a 2 Sin* } Y Cn) 


d* 


4 . 0 . 9 


4.0.4 


From 3. 0. 7 


du 


(n) 




- ' (•'Cot*}u (n, + {^)v ln) 

+ «‘«» * W ♦ <-EP > N** 0 ” 
(n) 


4 . 0. 10 


du 


Eliminating — between Equations 3. 0. 7 and 3. 0. 8, we obtain 


N_<n) . <<Azilf|ElCoy } u (n) ± { n( a 1 i i ^ )E '}v (n) - {^) E, }w (n) [. 


a 


♦ { v } N . 


X d * 

Eliminating [ — ^ 
we obtain £ 


1 1 

[ 3 between 3.0. 10 and 3.0. 11, 


M 


(n) B ^ (l-y 2 ) Dc °t^ j u (n) - ^ n(l-t/ 2 )P j y (n) 

1 a* a 2 Sin <j> 

+ t n 2 (l-u 2 )P } w (n) _ t (l- v 2 ) PCot<fr j x (n) 
a 2 Sin 2 <£ a 




Cn) 


+ { v } M 


Substitution of 4. 0. 11 and 4. 0. 12 into 4. 0. 9 yields 


Cn) 


dvj»„ v “ / .2 (n) 9 n / \ 'N 

• ~ { NSin <j> } 2_L_ + { n (l-v 2 )E , Cot0[ 1 +-^- ]}u W \ 
d$ 2 ” a 2 E' 


♦ [ 1+ a^r 1 > vln, 

+ { -nN + n (1 -v 2 )E?[I+ )„<">♦ { n <‘-^ - Pg?M } 

a 2 E' Sin 2 4 “ 

« 

+ { i/na} N .. W ? {im)M.. (n) -{Cot4}+ 0 (n) - {a 2 Sin* } Y W J 


4.0. 11 


4. 0. 12 


S 4.0.13 


4. 0. 5 


and substitution of 4. 0. 7 into 3. 0. 9 yields 


TTf--'- +l )u<n> +{Cot *> 


» 


+ ^E'(l-i/)Sin^ ^ ^ ^E'(l-V) * M i 


Cn) 


4. 0. 14 


Substitution of 3. 2. 12 and 4. 0. 1 into 4. 0. 14 leads to 


» 


^ j + D j dv v * w _ . f n r i _l D i \ ,.(n) 


TZ ' ST =t{ '3 iff? [lt ^, ))u 


a E 


+ {Cot^ [1 + -5- 1) v (n) T (^DCoti } w 
a^E' a E 1 Sin^ 

. f 2nD u (n) , i 2 \ ... W 

± { i[E T Br n* } + l E' i a-7)SiH? ) *0 


(n) 






J 


4. 0. 15 


Substituting 4. 0. 10 into 3. 0. 10, we get 


= {XHl— - ii + } w W - { i/Cot4< } 

aSin 2 * a 


_ / JL ) N < n > - (±) M tn) 

1 E' ' D ^ U 


Substituting 4. 0. 14 into 3. 0. 12, we obtain 






4. 0. 16 


M 


04> 


(n) , + ^ nD{l~i/)Co4 ^ w (n) 

2 0 . 


a Sin<£ (1+5 — ) 
a 2 E' 


4.0. 17 


- { nD(l-i/) 


aSin *(l+-5~^ 
a^E' 


}X (n) + { 


, ■ ■ ' 

a 2 E'Sin 4 (1+ ) 0 

a^E* 



Substitution of 4. 0. 4, 4. 0. 12 and 4. 0. 17 into 3. 0. 6 leads to 


—air 


<n) 


_ { U-y 2 )DCot 2 <j> ^ u (n) - ^ n(l-y 2 ) DCot^ ^ y (n> ^ 
" a Sin$ 


+ { n 2 D(l-,)Co^ j 2 + (Hv)] } w <n) 

a 2 Sin <£ (I+-5 — ) 

a^E' 

_ { liliLjP [ * 2n . l — + (l+y)Cot 2 0] )\ M 

a Sin <$(1+^9 — ) 
a 2 E' 

- {<!-./> Cot*} 


+ / 2 n D } . (n) 

“ a 2 E'Sin 2 *(l+-~ ) 0 

a^E' 


y 

( 4.0. 


18 


Upon substituting 4. 0. 17 into 4. 0. 4, we obtain 


(n>. { nfol- W . Co U. , n a D(t-> , ) , k W 

fi-f-2 — ) a 2 Sin 2 * (1+— 5 — ' 


a Sin *(l + ~9 — ) 


a E' 


nD 

a 41 E ' Sin 2 * ( 1+ 

a^E' 


- ^ 2 


)*« W + f^HT)n <n> 


J 4.0. 


19 




Substituting 4. 0. 14 into 4. 0. 7, we get 


~Y 


\\(n) 


N (n) a - { nD(l - v ) C ot * _ } w (n) { - nD (i:-^ } X 

a 3 Sin*(l+-y— ) a 2 Sin* (1+ /* ■— ) 

a E' aE' 


f 4.0.20 


+ ' aBi "*- U+ D_, 1 + ® ’ 


a^E’ 


We now 
taining 


BubBtitute 4. 0. U. 4.0.15. 4. 0. 19 and 4. 0. 20 Into 3. 0. 2, ob- 






2v Nii .>) 


) + T? 1)V 


f A. 0.21 


U-i/ 2 )E , Cot^ J . 2n DCotf tO-v>--Sil> w 

+ a a Sin ^ (1 +— -k ) 

a E* 


(n) 


< N+ -dr 


2 [(l-)- S,]> * W - «i-i')Cot*}N < 

Sin* ( I + —5 ) 


Cn) 


a E 


}+. <n> ±{-^V- t‘ + : 


2N 


+ ^ aSin <t> - a sin% 1 E* (1 — i' > (1+ — ^ 


u-ri } *0 


(n) 


a E* 


- ( a} X 


Cn) 


~J 


Substitution of 4. 0. 10. 4. 0. 11. 4. 0. 12. 4. 0. 16. and 4. 0. 17 into 4. 0. 6 
leads to 


4.0.8 


a i 2 

-J- » - {(l-i/)E'Cos*[(l+i/) + Si + n ^ D 1 ) u W 

E a 2 E'Sin 0 


+ { (1-t/ ) E'n[ (1 + i/) +~r + n ^ ~ y ~ 1 ) 

E a 2 E'Sin 2 0 


, r (l-i/ 2 )E'r o in 2. ,n 2 N . n 4 D 2n 2 DCot 2 0 


V 


4.0.22 


]}W 


to) 


r , , v . „ 4i r n 2 (l+ v) D A 2n 2 D 

{ (l-i/ )aE' Cot 0 [ — + 


a“E'Sin 4> a 2 E , Sin^^ 1 ^2~ ^ 


a ^ E ' ( 1 +— ^ — )(1 + i /) 
a E’ 

5 +|, Sin*l)X W 


a E* 


(n) N L^.ni 2 (n) 


- {(l+i/)aSin0}N. . + { a Su^+gj^ } M 


00 


+ ^ _2nDCot^___ j jj, to) _ { a? Sin 0 } Z* n) 

a 2 E' Sin 0(1+ -5 — ) 
a 2 E* 


J 


If we presume that D f E* and v are constants, upon differentiating 
Equation 4. 0. 15, we obtain 


9 


D , ,dv' n ' . r 2nD r Cos^+l_ n w (n) 


dX 


- f 4 n D Cot <f> ^(n) ^ r 2nD .... \ £2^ 
+ { aE’Sln 0 }X 1 1 aE'Sin* ' d* 


(n) 


d4* 


In) 


“ t + ^in«r * d 


e 


d<f> 


Now, substituting 4. 0. 10, 4. 0. 15 and 4. 0. 16 into 4. 0. 23, we get 




Cn) 


a E' d</» 

-{[-^ 4 -+i][i+V i} v<n) 

Sin0 a*E’ 


+ {*%££ i 1 + 


2 im 3 D 


^ 2nD f -,- 

Sin <T 1 J ‘ a 2 E'Sin 3 0 a E'Sin 3 * 


_ 2n(l + i/)D } ^(n) ? ^ 2nO(U i/)Cot _ } , 
a 2 E ' Sin ^ a E’Sin^ 


(n) 


, na . , D 2nD 

— ^ E ' Sin t ^ + a 2pV 


a^E’ a(E') Sin$ 


) " ‘ n, *<E H55* > M ( 


(n) 


+ {■ 


E'( 1- v ) Sint 


, 111 

% 


. 0. 23 


4. 0. 24 


Substitution of H«0»24 into k»0.13 yields 


4. 0. 10 




0 + 


2N 


di|i 


E’(l - 1/ ) ( 1 + ) 

a 2 E 1 


D ~dijr~ s - {nE'Cot^f 


, /, 2 W , . D V, ^ (n) , , n^E' r „ N 
+ ( 1 -„)( 1 + )j) u + <^ [fj, 

3L u 


+ (1- V 2 ) ( 1+ -B- » } v (n) + {nE'[v| 
a E' 


E 7 


♦ ^ - D + (i-, 2 )( i*sfo—a » > « (n) 

a^E'Sin* 


E ' a 2 E'Sin 2 <£ (1+5- ) 

a 2 E' 


y 


+ {nPOt^lCot [_2N + (l _„ ) J }x Cn) 

a E' (1+ 5—) 

a^E* 

+ (na [ v —• + pi s y-» J j N, , 

L E a 2 E‘(l+5- ) ** 

a^E’ 


± < n t ~7~T5 . }M ^ (n> * <Cot^}^ 0 (n) 


E’(l + -S ) 

a E' 

- { a 2 Sin<J> } 


J 


From 3. 0. 5, we have 


dM, 


(n) 


N < n) = - I Ui > - + » - m (n) - 

1N 0z a 3$ + aSin<£ m 00 a 


M, 


(n) 




V 


Upon differentiating 4. 0. 17, one gets 


4. 0. 25 


4. 0. 26 




dlV W _ -r / nD(i-v) (Cos 2 4>+1) ^ ... 

* 1 ~2i~srrr. dt i w 


(n) 


a Sin <t> (it — ) 
a^E' 

+ ^ 2n EKi-i; ) Coty j ^(n) 


aSin<£ ( 1 +-£— ) 
a^E' 


?{ _nD(l z£ ) , 

- o • _ i ✓ » . u \ 


dX 


(n) 


f 4.0. 


27 


aSin<£ (i +-g — ) 
a^E' 




- { ^Cot.^ p— >4>> 


+ { 


} 


dtj> 


(n) 


a^E'Sintf (1-h ^ — )’ T 6 
a^E* 


0 


a^E' Sin $(l+~ -P— -) ~ d * 
a^E* 


Substitution of 4. 0. 27 into 4. 0. 26 yields 


J 


N + { nD(l-y) (Cos^-fl) j w (n) 

a 3 Sin 3 * (1 + ) 

a^E' 




+ { 2 n D (1- v ) Cot* w Cn) + { nD(l-y) j 
1 2 r-, . _ , /, " D \ - _2 cij_ i / 1 . E J 


dX 


(n) 


a 6 Sin* (1 + A - -— ) 
a^E* 


+ { 


D Cot * 


a* Sin* (1 +~~ ) d* 
aE' 


s- rr- )* V M - <- 5— 2 D } 

a J E' Sin*( 1+ ) a d E'Sin* (1+^- ) 

m nl 


V 4.0 
(n) 


.28 


dvjj 


Q 


a E« 


a E* 


7 / n \ M M ( 2 Cot * ? 

+ ti^sr** M e& 1 / 


M, 


(n) 




4. 1. 1 


4. 1 Summary of First Ord er Equations 

From the previous section we summarize the results as follows: 


- (*«*) « to) T v<">. «Hrt) 


♦ t* > N W (n> 


r 4. 1. 1 


dv^_ t f n \ ..(n) A \ ..GO 

" ± { -snrr } u + {Cot * } v 


+ { 2nDCoy _ } w (n) +{ 2nD_ 


a^E' Sin^ (1-t — 4r- — ) 
a^E' 


a E' Sin $ (1 -h-s ) 

a Z E’ 


B v )X<n) ^ 4. 1.2 


+ { 


XT v }4 o 


(n) 


E'(l-v) Sin <f> (1+^^-) 
a E’ 


$ w -* w 


4. 1. 3 


2 

l 

c 

a Sin 4 * 


. {™- - ( i±i.)}w fa> - {^Cot*)X (n) 
d 4> - 


- <B*> N # <n) - <IT> M ,j 


(n) 






4. 1. 2 


. ^[(.-^cotV n a-4 -U)« w 

Sin <f> 


i ‘-TigS 4 - i “-*-*>♦ > v<n> 


+ ( -( ‘ _ -J' i 2 > . e ' Cot* , 2n 2 DC0t [(1 . y) >) 

a a ain ^(1 +J-) E > 415 

2n 2 D[(l- I /)- g,] a E ' 

~ {Nt °V X 

a 2 E' 


* <^> +* <n> ±<-V t »♦— J 

v v aSin $ k • n — i 


a^E’ 


- {a} X 


(IM O n 

-W o- i U-t'JDCo \ „(n) — s n(l-v 2 )DCot<i> ' (n) 

do 1 5 -f u + t — / v 

a ^ Sin 0 

j _ / a D(l — i/) Cot $ r tt; /i v -i ■> (n ) 

— 1 <»♦-&•> + < ’ 1)w 
a E' 


- < ^ I z^fr -d- ♦ 0 - *> cot 2 * ] } *<"> 

Sin 0(1 + — w — . ) 

a^E' 

i (n) 

-«l-rtC**)M ^ +{s5r ) ♦ 


+ { — o 


2 n D 


a E' Sin* 


-5-/ V 


a E' 


4. 1. 6 


d4> 


(n) 


^rJ- = -{(1-1/)E' Cos 4> [ (H-W + g, + n . ] } u Cn> 

d< P E a 2 E' Sin $ 


+ { (1-y ) E'n[ (1 + ^N + n(l ±i ^ ] } y (n) 

E a Z E‘Sin<* 


(1-i/ 2 )E« rcA, n 2 N 


n 4 D 


, / iej r ci>j . 11 a” *. “ u 

+ < -sir ~ [Sin ^ETrnr + a -2iTir2^ 


. 2 n 2 D Cot 2 <b , „.(n) 

+ — n g — / w 

a 2 E' (1+v ) (1+ -#-) 


9 ‘ 
a Z E' 


n 2 (i w ) D 2 n 2 D 


a E* 


- {(1- 1 / ) aE' Cot$ [ + 

a E Sin$ a^E' Sin<£0 +- o ) 

+ £ I Sini]}X (n) - {(l+i/)aSin*> N^^ <n) 


t { g- a 2 Sin* + ^ } M. + ^aDCoy } 4 , 

D SIS? M ~ a^E'Sin? (l+-£— ) 

a 2 E' 


- { a 2 Sin <j> } Z ^ 


d4< 
d <j> 


e 


= i 


1 + 


2N 


r v N 

+ {nE'Cot 4> [ -fjT- 


E'(l-i') U + "F 
a^E' 


2 • D o i (n) . n 2 E 1 r vN 

+ ( l-v 2 ) ( 1+ —rr~ )] ) u Sin<t> 1 E' 

a 2 E’ 


+ (1- v 2 ) (1 + -J^-) ] ) v(n) + (nE* 

a 2 E' 


2 

2N P(l+yn 2 -(l+i/) Sinf) + (i- v 2 ) (1 + P --)] ) 

^ a^E'Sini (1+ -S-) a 2 E' Sin ^ 

a*E' 


f nD(l-M/)Coti r 2N + (l-i,) ] } \ <n) 

- 1 a E'(i+-#-) 

a 2 E' 


+ { na [ v 


N . 2N 
E 1 E' 


D 


a ^E'O* -j— > 
a 2 E* 



- t n[ ^ + »1}>« W W - {Cot* ) +„ (n) 

e'(i+_JB ) 

a 2 £’ 

- { a 2 Sin <f> } Y (n) U 



4. 1. 5 


<„> . ( w- 2 )g. s at.) u w ± iidfegp , v<"> 


+ { v) N 


(n) _ { ( l-v)E ' } w (n) 


(n) . _ {< iz y|DCot^ ^ u (n) ^ { n(l-^ 2 )P } y (n) 


a 2 Sin 4 


+ { n 2 (l-»/ 2 )D } w (n) _ { (1- i/* ) D Cot <f> } ^(n) 

a 2 Sin 2 * a 

+ M M < n) 


N ("> * + { ."Pg r^ cg^L > w (n) ± { 


a°Sin0(l+ “ 5 =- — ) 
a^E' 


a * Sin</> (1 +— s — ) 
a^E' 


1 ii ' • « 

a Sin0 { 1 + — s — ) 
a E 1 

Mfl W. +{ »DtjiCoy^_ )>> + ( nD( ;- -ii> } x< 

a 2 E' Sin0(l+ ) aSin<^ (1 +-* — ) 

a 2 E* a E' 


T \ -p f) ' 

a^E' Sin* (1+ — j- ) 

a^E' 

n a W . { ) w W - { } * (n) 

a 3 Sin 2 0 ( 1 +-£■—) a 2 Sin 2 0(l+-o — ) 

a 2 E' a 2 E' 


} a‘SIA 0 } *4 


a 3 E* Sin 2 * ( l " ® 1 aSIn* 0 


a*E’ 



4. 2. 1 


* 
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< 


* 
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4. 2 Nondi mensionalizat ion of Equations 


For engineering problems, it is frequently advantageous to present 
results in a nondimensional form. This is particularly true where design 
techniques or general results are involved. The nondimensional form per- 
mits one set of data or calculations to be applied to a large number of prob- 
lems. Examination of the equations 4. 1. 1 through 4. 1. 14 shows that the 
selection of certain quantities as the nondime nsionalizing factors will re- 
sult in some simplification of the equations. Accordingly, equations 4. 1. 1 
through 4. 1. 14 have been nondimensionalized. These nondimensional equa- 
tions are presented in Chapter 1 as equations 1. 1. 1 through 1. 1. 16. 

The scheme of nondimensionalization is presented below. In general 
"a primed quantity" will denote a nondimensional quantity. There is one ex- 
ception. This is E', a dimensional quantity as defined previously. 


(n>« 

u 

u (n) 

3 _ 

v <n)’ = 

>> 

a 


a 

(n)’ 

w 


** 

If 


a 

E' 


. < n) 


H (n) 

- ~sw 

aE» 

V n) ' 

*0 

aE r 

X< n) ' = 

a „<n) 
E' A 

Y <n)' 

a v (n) 
- eT 1 

z (n)‘ - 

a ~(n) 

w z 

N' 

N 

= W 

D* 

D 

a 2 E' 


"N 


> 4. 2. 1 




5.0. 1 



Integration of Equations 

The solution of the problem for stresses and deflections in a 
spherical shell has now been reduced to the integration of a set of eight 
linear first-order differential equations. This set is readily integrated 
on a digital computer using one of several integration procedures. The 
program prepared as part of this project utilizes a fourth-order Runge- 
Kutta process. 

The numerical integration presents several mathematical compli- 
cations, each of which has been resolved in the accompanying program. 

The first, and most easily overcome, complication results from 
the fact that the problems to be solved are boundary value problems where- 
as numerical integration schemes are directly applicable only to initial 
value problems. Since the governing differential equations are linear, 
the boundary value problems can be solved by constructing a linear com- 
bination of solutions to judiciously chosen initial value problems so that 
both the initial and final boundary conditions are met. 

The second difficulty occurs when the angles ^ = 0 and/or $ = jr 

lie within the region of integration. The differential equations contain 

singularities at these points. In order to avoid these singularities, such 

artifices as a small hole, a small rigid plug, or a small elastic plate or 

cap may be introduced into the program. In the case of the program pres- 

-6 

ented, an extremely small hole subtending a half-angle equal to 10 times 
the outside angle is used. 

The program may be readily modified to use any of the alternate 
possibilities which have been mentioned above. 


